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1 = 3/2 and 7rA* — > pA for 1 = 3/2. As the temperature increases, the rise or 
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loosened bound states of final mesons, and the total-mass difference of the initial 
and final mesons. The temperature-dependent cross sections and meson masses are 
parametrized. 
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1. Introduction 

In spite of scarce experimental data of cross sections for inelastic meson-meson scat- 
tering, these cross sections can be calculated from quark potential models, parton dis- 
tributions or effective meson Lagrangians that respect various symmetries. This has 
been seen from the dissociation cross sections of J/ip in collisions with mesons, which 
have received much attention. For example, cross sections for reactions like n + J/ip — > 
DD* + D*D + D*D* and p + J/ip -> DD* + D*D + D*D* were obtained with the 
quark-interchange mechanism [1,2] in different quark potential models [3-6]. The quark- 
interchange models give the characteristic that cross sections for endothermic reactions 
first increase from threshold energies and then decrease as the center-of-mass energy of 
J/ip and hadron increases. A very small low-energy nucleon- <//■?/> break-up cross section 
was obtained by using the operator product expansion for the elastic scattering amplitude 
of J/ip [7-9]. The dissociation cross section for n — J/ip : evaluated in QCD sum rules [10] 
in the soft-pion limit, increases gradually in contrast to rapid growth near the threshold 
energy obtained in meson exchange [11-15] and quark interchange models [3-6]. 

In comparison to meson- J/ip reactions, however, cross sections for inelastic scattering 
of a light meson by another light meson are rarely studied. Inelastic scattering between 
light mesons does occur in hadronic matter and the cross sections involved influence the 
time dependence of meson momentum distributions and flavor dependence of the measured 
momentum distributions at kinetic freeze-out. In order to understand this influence, the 
cross sections for inelastic scattering and their relevant characteristics must be studied. 
Since hadronic matter produced in Au-Au collisions at the Relativistic Heavy Ion Collider 
mainly consists of pions, rhos and kaons [16-23] , in this work we pay attention only to the 
nonresonant reactions of the four mesons n, p, K and K*, which are taken to be governed 
by quark-interchange processes. In Ref. [24] we have calculated in-vacuum cross sections 
for the nonresonant reactions — > pp for 1 = 2, KK — > K*K* for I — 1, KK* — > K*K* 
for I = 1, nK -> pK* for / = 3/2, irK* -> pK* for / = 3/2, pK -> pK* for / = 3/2 and 
nK* — > pK for I = 3/2. The cross sections for the seven endothermic reactions depend 
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on the center-of-mass energy y/s of the two initial mesons and the energy where the 
maximum of cross section occurs is mainly determined by the maximum of | P' \ / s \ P |, 
where P and P' are the momenta of the initial and final mesons in the center-of-mass 
frame, respectively. The endothermic reactions have maximum cross sections ranging 
from 0.47 mb to 1.41 mb. These cross sections are obtained with a quark-quark potential 
that includes the linear confinement and one-gluon-exchange potential plus perturbative 
one- and two- loop corrections [6]. At nonzero temperature the linear confinement is 
modified to become weaker. At a temperature below the critical temperature T c of the 
QCD phase transition, medium effects show up in the region where the quark- ant iquark 
distance is larger than 0.3 fm [25]. When the quark-antiquark distance is large enough, the 
quark- ant iquark potential at a given temperature becomes constant. This temperature- 
dependent potential must change the wave function of quark-antiquark relative motion in a 
meson and cross sections for inelastic meson-meson scattering are expected to change with 
temperature as well. Therefore, in this work we study the dependence of cross sections 
on the temperature for the nonresonant reactions nn — > pp for 1 = 2, KK — > K*K* for 
1=1, KK* -> K*K* for I = 1, nK -> pK* for / = 3/2, ttK* -> pK* for I = 3/2, 
pK — > pK* for / = 3/2 and nK* — > pK for / = 3/2. Until now the temperature 
dependence of the cross sections has not been studied in experiments or theory. 

In the next section, we present a potential of which the short- distance part is given by 
perturbative QCD and the large-distance part is displayed in the lattice gauge results of 
Ref . [25] . The Schrodinger equation with the potential is solved to get the quark-antiquark 
relative-motion wave function of a meson. A convenient framework for application of the 
temperature-dependent potential to inelastic meson-meson scattering is the Born approx- 
imation for the quark-interchange processes. In Section 3 we review formulas of cross 
sections for meson-meson nonresonant reactions that are based on the Born approxima- 
tion. In Section 4 we show numerical results of cross sections and give relevant discussions. 
Parametrizations of the cross sections are given. Conclusions are in the last section. 

2. Potential and wave functions of quark-antiquark relative motion 
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In Ref. [25] the heavy quark potential was assumed to be equal to the free energy of 
heavy quark-antiquark pair and the lattice calculations provided the temperature depen- 
dence of the potential. When the distance r between quark and antiquark is large enough, 
the potential at a given temperature exhibits a constant value and becomes a plateau. 
This value depends on the temperature T and decreases with increasing temperature. 
The constant value at large distances can be parametrized as 



V ab (r) = -— ■ --D 



T x 4 



(1) 



with D = 0.7 GeV and T c = 0.175 GeV. Here A are the Gell-Mann matrices for the 
color generators. D is a fit parameter, but happens to equal the height of the plateau of 
T = 0.74T C . 

The relative color orientation of two constituents (quarks and antiquarks here) can 
be indicated by wave functions of the two constituents. The expectation value of 4f • 4r 
depends on the relative color orientation, so does the heavy quark potential given in Eq. 
(1). The dependence of the potential on the relative color orientation is supported by 
hadronic physics and lattice calculations. We discuss evidence in the three cases: (a) 
T = 0, (b) T > T c and (c) < T < T c . 

(a) Since the establishment of QCD the dependence of confinement with 4f • y on the 
relative color orientation has been widely used in hadronic physics, for instance, meson- 
baryon scattering [26] and meson-meson scattering [5]. The expectation value of 4f • 4f 
is negative in the color singlet state and positive in any color octet state. This complies 
with the fact that the color singlet state is observed and the color octet states are not 
observed. 

(b) When temperature is above T c , it is expected in Ref. [27] that the interaction in 
an antitriplet diquark system is only half as strong as in a color-singlet quark-antiquark 
system, i.e. obeys Casimir scaling, and the interaction depends on the relative color 
orientation. 

(c) When temperature is between and T c , the situation is complicated. From Figs. 
15, 17 and 20 in Ref. [28], we can derive that the color-singlet free energy does not equal 



the color-octet free energy at large distances in the region 0.907T C <T <T C . It can also 
be expected from Fig. 3 of Ref. [27] that the interaction in an antitriplet diquark system 
is nearly half as strong as in a color-singlet quark-antiquark system, i.e. the interaction 
depends on the relative color orientation in the region 0.87T C <T <T C . 

Much attention has been paid to construction of the color octet potential from per- 
turbative QCD and lattice data in the study of hybrid mesons. The color octet potential 
relies on the spatial extension of gluon field confined in hybrid mesons. In the string 
description the symmetry of the gluon field is labeled by the spins £, II and A about the 
axis connecting the quark and antiquark, a PC value and an additional reflection sym- 
metry for E states. It was obtained in vacuum that the color octet potential is different 
with respect to different symmetries of the gluon field [29,30]. Temperature dependence 
of the gluon field affects temperature dependence of the color octet potential. How does 
the gluon field confined in a meson depend on temperature at < T < T c ? This is a 
difficult problem. 

In summary, we still lack of the lattice study of the color octet potential at < T < T c 
and the heavy quark potential depends on the relative color orientation. 

The large-distance plateau of the heavy quark potential indicates the onset of string 
breaking at a certain distance [27]. The string breaking occurs when the energy stored in 
the string exceeds the mass of a light quark-antiquark pair [31]. The combination of the 
light quark and the heavy antiquark and the combination of the light antiquark and the 
heavy quark form open heavy flavors. When the temperature is higher, light mesons in 
medium more effectively take a kind of flip-flop recoupling of quark constituents [32-34] 
and string breaking occurs at a shorter distance. Hence, the plateau appears at a shorter 
distance and confinement gets weaker. Consequently, the wave function of the heavy 
quark-antiquark pair becomes wider in space. 

However, uncertainty of the heavy quark potential exists because the potential is de- 
rived from Polyakov loop correlation functions. The Polyakov loop at the space coordinate 
fi is indicated by L(fi). The Polyakov loop correlation functions are related to the free 
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energy of heavy quark-antiquark pair F(T, r) 



Tln< L(fi)L+(f 2 ) >= F{T,r) + C 



(2) 



where r =| fi — r 2 | and C is a normalization constant. The free energy leads to the 
internal energy by 



The entropy S(T,r) = —dF(T,r)/dT is independent of r at large distances and depends 
on r at the other distances. If TS is small, U(T,r) ~ F(T,r); otherwise, the internal 
energy deviates from the free energy. Taking the free energy as the heavy quark potential 
like Ref. [25] is an approximation. Then the internal energy was suggested in Refs. [35,36] 
as the heavy quark potential. Very recently the heavy quark potential is argued to be the 
internal energy which is the expectation value of the difference of the Hamiltonians with 
and without the heavy quark-antiquark pair at rest [37,38], and this results in dissociation 
temperatures of heavy quarkonia that agree quite well with the values from lattice studies. 

How to find the correct heavy quark potential is an important problem but has not 
been solved so far. It has been proposed that the heavy quark potential takes the form 
£U(T,r) + (1 — £)F(T,r) where the quantity £ is between and 1 and was determined 
in models [39-45]. Since the present lattice calculations provide the internal energy U 
that includes the internal energy of the heavy quark-antiquark pair and the gluon in- 
ternal energy difference U g (T,r) — U g0 (T) where U g (T,r) and U g0 (T) correspond to the 
gluon internal energies in the presence and absence of the heavy quark-antiquark pair, 
respectively, it was proposed by Wong [41] that the heavy quark potential should be 
U(T,r) — [U g (T,r) — U g o(T)]. In the local energy- density approximation that adopts an 
equation of state for the quark-gluon plasma, the heavy quark potential can be repre- 
sented as a linear combination of F and U . With F and U obtained in lattice calculations 
in quenched QCD, the potential gives spontaneous dissociation temperatures of heavy 
quarkonia that agree with those obtained from spectral analyses in quenched QCD [43]. 
This is an improvement in comparison to the use of F as a potential. The work of Wong 
implies that the heavy quark potential properly defined gives reliable results. Even though 



U(T,r) = F(T, r) + TS(T, r) 



(3) 
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the method of Wong has not been applied to hadronic matter, we still expect that the 
heavy quark potential for hadronic matter, if obtained, can produce more reliable cross 
sections for the meson-meson nonresonant reactions in the nonzero temperature region of 
hadronic matter than those from F. 

At very small distances r < 0.01 fm, the potential obtained from one-gluon exchange 
plus one- and two-loop corrections in perturbative QCD is [46] 

462 In w" 



53 \ 1 
2-Ye + ■==)- 



(4) 



75 J w 625 w 

where 7^ is the Euler's constant and w = \n(l/A^r 2 ) with the QCD scale parameter 
Aj^jg determined by Eq. (2.13) in Ref. [46]. 

An interpolation between the constant confinement at large distances given in Eq. 
(1) and the spin-independent perturbative potential given in Eq. (4) produces a central 
spin-independent potential 



V ab (r) = -— ■ --D 



ta„h(Ar) + ^4|^>exp(-Br), (5) 



where 



and 



A = 1.5[0.75 + 0.25(^) 10 ] 6 GeV, 



are fit parameters, and 



E = 0.6 GeV, 



A = v^&oASTrV, 



(6) 



(7) 



(8) 



a' = 1.04 GeV" is the Regge slope and b = 11 — |jV/ with the quark flavor number 
Nf = 4 [46]. This potential is different from the parametrizations given by Digal et al. 
and Wong [47] . The potential contains the dimensionless function [46] 

Ab [°°dQ ^ 2 K Q 
^ Jo Q Q 2 A 



(9) 



with K = 3/16rr 2 a', where Q is the absolute value of gluon momentum Q and p(Q 2 ) is 
given by Buchmuller and Tye [46]. The quantity p—-^k arises from one-gluon exchange and 



perturbative one- and two-loop corrections. exp(— Er) is a medium modification factor 
to the potential of one-gluon exchange plus perturbative one- and two-loop corrections. 
The temperature correction to the one-gluon-exchange potential with the limit shown 
in Eq. (4) is the difference between the second term in Eq. (5) and the perturbative 
potential 4f • ^ll [46] ■ The temperature correction is completely negligible at very 
short distances and obvious at intermediate and large distances. 

The parametrization in Eq. (5) versus lattice gauge results is plotted in Fig. f . It is 
clearly seen that a lower plateau at large distances corresponds to a higher temperature. 
Plateaus at T/T c = 0.97, 0.94, 0.9, 0.84 approximately begin at r = 1. 15 fm, 1.18 fm, 1.26 
fm, 1.38 fm, respectively. Hence, a higher plateau begins at a larger distance. Confinement 
can be assumed to be flavor-independent in hadronic physics. For example, it was shown in 
Refs. [48-53] that quark-quark potentials with flavor-independent confinement, a Coulomb 
term and hyperfine interactions can consistently describe a large body of data like masses 
from light to heavy hadrons. The light hadrons may consist of only up and down quarks 
and the heavy hadrons may contain charm and bottom quarks. The flavor dependence of 
the quark-quark potentials is relevant to quark masses in the hyperfine interactions. The 
success of the potentials renders that the flavor independence of confinement is universal 
and the hyperfine interactions must be flavor-dependent. Therefore, the potential in Eq. 
(5) is reasonably flavor-independent. The potential obtained by Karsch et al. for heavy 
quarks is applied to light quarks. It is shown by the lattice calculations that screening sets 
in at distances r 0.3 fm. However, present lattice calculations at finite temperatures 
are probably not yet precise enough to reach sufficiently short distances. Hence, a lattice- 
based potential has a degree of freedom in choosing its form at short distances. We 
construct the potential in Eq. (5) from the fact that the quark-quark interaction at 
r < 0.01 fm originates from one-gluon exchange plus loop corrections in perturbative 
QCD. The degree of freedom is removed at r < 0.01 fm. 

The medium screening obtained in the lattice gauge calculations at present affects 
only the central spin-independent potential. This allows us to keep using the spin-spin 
interaction relevant to perturbative QCD, as done by Wong [47]. In our work, the spin- 



8 



spin interaction arises not only from one-gluon exchange but also from perturbative one- 
and two-loop corrections [6] 

A a \ b 16n 2 3 s a ■ s b \ a \ b An 1 d 2 v(\r) s a ■ s b 

V ss = -— ■ ———5%r) + — — — , 10 

2 2 25 m a m b 2 2 25 r dr z m a m b 

where s a and m a are the spin and mass of constituent quarks or antiquarks labeled as a, 
respectively. This expression of the spin-spin interaction comes from the application of Eq. 
(7e) in the transformed Hamiltonian obtained by Chraplyvy [54] from the two-constituent 
Hamiltonian that includes the relativistic potential originating from one-gluon exchange 
plus perturbative one- and two- loop corrections [6]. The transformed Hamiltonian was 
obtained from an application of the Foldy-Wouthuysen canonical transformation to a 
relativistic two-particle Hamiltonian. The spin-spin interaction is related to the terms of 
the direct product of two Dirac a matrices [55] in the relativistic potential. 

Given the masses m u = = 0.32 GeV for the up and down quarks, the Schrodinger 
equation with the central spin-independent potential in Eq. (5) is solved at T = to 
obtain a radial wave function R q q(T = 0,r) for the quark- ant iquark relative motion of 
7r and p mesons. Assuming all the mesons in the ground-state pseudoscalar octet and 
the ground-state vector nonet taking the same spatial wave function of quark- ant iquark 
relative motion as the n and p mesons, the spin-spin interaction leads to the mass splitting 
between a pseudoscalar meson and a vector meson with the same isospin quantum number 

< v - >= r£k^ iT = »■ r = 01 - f *' r ^*«< T - <»> 

where R q q satisfies the normalization condition J °° drr 2 R 2 g(T,r) = 1. At m s = 0.5 GeV 
for the strange-quark mass, the spin-spin interaction yields the mass splittings m p — m n = 
0.5989 GeV, m K * — m K = 0.3833 GeV and |m w + \m < $ > — m v = 0.3622 GeV, where rrii 
(i = 7r, p, K, K*, T), ou, 0) represent the masses of ir, p, K, K*, r], uj and 0, respectively. These 
mass splittings can be compared to the experimental values 0.6304 GeV, 0.3963 GeV and 
0.3930 GeV, respectively. 

The spin- averaged mass of a spin-0 meson and a spin-1 meson with the same isospin 
is one-fourth of the spin-0 meson mass plus three- fourths of the spin-1 meson mass [3]. In 
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Table 1 we list vacuum masses of 7r, p, K and K*, the spin-averaged mass m np of ir and 
p and the one Wikk* of K and K* . The theoretical values of m np , Wlkk*-, m P , m K and 
rriK* approach the corresponding experimental data. The pion mass from our calculations 
almost doubles the experimental value. These are understandable. The potential given 
in Eq. (5) has the behavior of tanh(Ar) at large distances and cannot mimic the linear 
confinement. In vacuum the pion with the lightest mass among mesons has the largest 
radius and is sensitive to the potential behavior at large distances. Therefore, the exper- 
imental datum of pion mass must differ from the value derived from the potential. For p, 
K and K* with masses larger than 7r, less sensitivity to the potential behavior at large 
distances does not lead to a large separation of the experimental and theoretical masses. 
If the theoretical pion mass is used, cross sections for pions in collisions with mesons 
are not reliable. Hence, the experimental values of meson masses are used in calculating 
meson-meson cross sections at T = GeV in this work and, as will see, the cross sections 
are similar to those obtained from the potential with the linear confinement in Ref. [24]. 

If the quark masses are set equal to zero, the QCD Hamiltonian is symmetric under the 
chiral group SU (3) x SU (3). Although the spontaneous breakdown of the chiral symmetry 
is widely expected, the problem of calculating the spontaneous chiral symmetry breaking 
in QCD has not yet been satisfactorily solved [56,57]. However, the spontaneous symmetry 
breaking can be related to the quark condensate which vanishes in perturbative QCD 
and does not equal zero in the non-perturbative region. With the explicit breakdown 
of the chiral symmetry due to the current quark masses, the square of pion mass is 
proportional to the product of the quark condensate and the sum of the up and down 
quark masses at the lowest order in chiral perturbation theory. Therefore, the small pion 
mass comes from the small quark masses and the non- vanishing quark condensate in the 
non-perturbative region where confinement sets in. While only the spontaneous symmetry 
breaking occurs, Goldstone bosons and quarks are massless and the soft-pion theorems 
obtained in the vanishing pion four- momentum are exact [58]. While the explicit chiral 
symmetry breaking also occurs, bosons get masses and the soft-pion theorems become 
approximate. The deviation of the theorems from the experimental data has been studied 
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in chiral perturbation theory and the soft-pion theorems have been corrected [59-62]. In 
our work all the mesons in the ground-state pseudoscalar octet and the ground-state 
vector nonet are assumed to take the same spatial wave function of quark- ant iquark 
relative motion as the pion and rho mesons, but the quark masses are the constituent 
masses which are much larger than the current quark masses used in chiral perturbation 
theory. The constituent quark masses show the explicit chiral symmetry breaking so that 
mesons get masses and the soft-pion theorems are violated. The cross sections that we 
will obtain from the constituent quark masses and the nonzero meson masses must deviate 
from the cross sections obtained from vanishing masses of quarks and mesons. Therefore, 
the deviation of our results from the soft-pion theorems should be expected. 

3. Formulas for cross sections 

Let TOj and P« = (Ei, Pi) be the mass and the four-momentum of meson i (i — 
Q1Q1, Q2Q2, Q1Q2, Q2Q1), respectively. The Mandelstam variables for the reaction A(qiq 1 ) + 
B(q 2 q 2 ) -»• C(qiq 2 ) + D(q 2 q 1 ) are s = (E qiqi + E q2q2 f - (P qiqi + P q2q2 f and t = (E qiqi - 
Eqiq2) 2 ~(Pqiqi ~ Pqifc) 2 ■ I n the center-of-mass frame the meson A{q{q l ) has the momentum 

— » — * — * — * 

P = P qiqi = —Pq 2 q 2 an d the meson C(q{q 2 ) has the momentum P = P qiq2 = —Pq 2 q 1 - P 

— * 

and P' are expressed in terms of s by 

l^(v^)| 2 = ^ { [s ~ (m 2 qiqi + ml q2 )] 2 - Aml- qi ml q2 } , (12) 

l^(v^)| 2 = I { [s ~ K, 2 + <J] 2 - 4< ? - 2 <4 . (13) 

Denote the angle between P and P' by 6. The cross section for A(qiq 1 ) + B(q 2 q 2 ) — > 
C(q{q 2 ) + D(q 2 q l ) is 

^= qo 1 7^7^ r de\M & (s,t)\ 2 sine. (14) 



32tts \P(^7s)\ Jo 

This formula provides the dependence on the total energy ^/s of the two initial mesons 
in the center-of-mass frame and is valid for the interchange of the two quarks (qi and 
q 2 ) or of the two antiquarks {q~i and q 2 ). The interchange of quarks brings about two 
forms of scattering that may lead to different values of the transition amplitude The 
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forms are known as the prior form and the post form [63-65] . Scattering in the prior form 
means that gluon exchange takes place prior to the quark or antiquark interchange. The 
transition amplitude in the prior form is [24] 



Ayi prior _ in 771 <Tp <Tp 9 771 — f d 3 p qi g 2 d 3 p 
JVl fi - V ^qiqi^qril^qi&^qiqi J (27f) 3 



' < Pqiq2^Pqiq2) 1 Pq2qi^P>12qi)^Yqiq2 ^9192 "I" ^9192 ^9192)09191 (^*9l 9l) ^92 92 (.P9292)' 

(15) 

where ipabiPab) is the product of color, spin, flavor and momentum-space wave functions of 
the relative motion of constituents a and b and satisfies / ^^^^ipab^abipab) — -• The 
relative momentum of a and 6 is p a b- Scattering in the post form means that the quark 
or antiquark interchange is followed by gluon exchange. The transition amplitude in the 
post form is [24] 



' 091 92 (Pll 92 ) 0V; 91 (^92 91 ) ^91 91 091 91 (^91 91 ) 092 92 (^92 92 ) 



^Pqiqi d3 P 9192 , 



(2tt) 3 (2tt) 3 



^ j "(271^)^ (27Ty^ _ ^'^ L ^ 2 52 ) V^q^ (Pg25l )^ / 52<?2^/ ; giQl (Pgi^l )V'l32Q2 (^?2Q2) 

+ j ^2tt) 3 (27r) 3 ^ qi ^ 2 92 ^ ^9- 9 1 91 ) ^91 92 091 91 (^91 91 ) 092 92 (^92 92 ) 

^ J ~(27r)^ (l^Try^ - ^^^ 2 52 ) V'q^ gi (^?25i )^ig2^giQi (^<?igi )^g2Q2 (^2^2) 



(16) 



The transition amplitudes in the prior form and in the post form are equal to one 
another when the potential and the wave function of quark-antiquark relative motion are 
those used in the Schrodinger equation [63-65]. Otherwise, M^ lor 7^ M^ st for inelastic 
scattering. The inequality yields different cross sections corresponding to the two forms, 
which is the so-called post-prior discrepancy [63-65]. 

In the Schrodinger equation we have only used the central spin-independent potential 
in Eq. (5) since the spin-spin interaction in Eq. (10) contains the delta function that 
can not be correctly dealt with in the equation. But in the transition amplitudes where 
the Fourier transform of the spin-spin interaction can be correctly dealt with, we use 
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the Fourier transform of both the central spin-independent potential and the spin-spin 
interaction: 



Kb(Q) = -y-yi D 



1.3- 



(2vr) 3 5 3 (g) - ^ / dr- 



exp(2Ar) + 1 



2 2 7 * (£ 2 + Q 2 + g 2 ) 2 - 4QV 

A a A 6 16tt 2 s a • s fe A a A 6 16?r 2 A f°° d 2 v (x) . / Q \ s a • s 6 
+ "Tr • — — ^ / ax — — — sin —x 



2 2 25 m a m(, 2 2 25Q 7 rfx 2 \ A / m a mb 

(17) 

Therefore, the post-prior discrepancy occurs in our calculations. We take the average 
of the cross section in the prior form and the one in the post form. Each of the two 
cross sections related to the prior form and the post form is the unpolarized cross section 
obtained from the cross section in Eq. (14) 

CTU " POl( ^ = (2 S , + l)W B + D ^ (2S+1) ' T(S ' m " V5) - <18) 
where Sa and S B are the spins of A and B, respectively, and S is the total spin of the 
two mesons allowed by the reaction A + B — > C + D. The cross section o- un P o1 (-^/s) is 
independent of the magnetic projection quantum number ms of S. 

4. Numerical results and discussions 

The quark masses m u = m A = 0.32 GeV and m s = 0.5 GeV determined in the fit to the 
experimental mass splittings in Section 2 are kept unchanged in hadronic matter. They 
are used in the Schrodinger equation with the central spin-independent but temperature- 
dependent potential in Eq. (5). The lowest-energy S-wave solution of the Schrodinger 
equation is a temperature-dependent radial wave function of the quark-antiquark rela- 
tive motion of mesons in the ground-state pseudoscalar octet and the ground-state vector 
nonet. When the temperature increases, the peak of r times the .S-wave radial wave func- 
tion, rR q q(T,r), moves to larger quark-antiquark distances and the meson's root-mean- 
square radius increases. This reflects the phenomenon that with increasing temperature 
any bound state becomes looser and looser while confinement gets weaker, i.e. the poten- 
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tial plateau at large distances decreases with increasing temperature and a higher plateau 
begins at a larger distance. 

Another consequence of the temperature-dependent potential in Eq. (5) is that masses 
of 7r, p, K and K* decrease with increasing temperature in the region 0.6 < T/T c < 0.99. 
The mass splitting of a spin-0 meson and a spin-1 meson with the same isospin is calculated 
according to 

< V - >= 75^1^ '' = °> - f *^^,r)]. (19) 

The spin-averaged mass of a spin-0 meson and a spin-1 meson with the same isospin 
is one-fourth of the spin-0 meson mass plus three-fourths of the spin-1 meson mass [3]. 
The spin-averaged mass of n and p equals the sum of quark mass, antiquark mass and 
the energy of the relative motion obtained from the Schrodinger equation. Since for the 
quark- antiquark relative motion we take the same spatial wave functions of K and K* as 
that of 7r and p mesons, the spin-averaged mass of K and K* equals the sum of quark 
mass, antiquark mass and the nonrelativistic-Hamiltonian expectation value of the wave 
function. After mass splittings and spin-averaged masses are obtained, we find meson 
masses of which the mass of the spin-0 meson is the spin-averaged mass minus three- 
fourths of the mass splitting of the spin-0 and spin-1 mesons, and the mass of the spin-1 
meson is the spin-averaged mass plus one-fourth of the mass splitting. The temperature- 
dependent meson masses are plotted in Fig. 2. The reason for the falloff of masses with 
increasing temperature is that the first term of the potential related to the large-distance 
plateau gives a smaller contribution at higher temperature. When T — > T c , the masses of 
the 7i and p mesons approach and 0.006 GeV, respectively. This indicates that the n 
and p mesons are almost massless at a temperature very close to T c . In addition to the 
falloff of masses, the mass splitting of n and p and the one of K and K* are both close to 
zero at T — > T c . When T — > T c , the masses of K and K* approach finite values 0.179 GeV 
and 0.183 GeV, respectively. The kaons and vector kaons become degenerate in mass near 
T c . Similar to Ref. [3], a constant of 0.88857 GeV is subtracted from the energy of the 
Schrodinger equation for n and p and from the nonrelativistic-Hamiltonian expectation 
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value of the wave function for K and K*. This subtraction makes the theoretical values 
of the spin-averaged masses at T = approach the experimental data and the masses of 
it and p go to zero at T — > T c . The subtraction does not influence the sizes of meson 
bound states and the mass splittings. 

The meson masses in units of GeV are parametrized as 



m„ = 0.41 



1 - 



T 



1.05T,. 



11.88 



3.81 



m p = 0.7 



ttik = 0.63 



' - 1 T, 



4.29 



1.14 



1 - 



T 



1.157; 



9.31 



4.4 



ttik* 



0.84 



1 - 



T 



1.057; 



4.16 



(20) 

(21) 
(22) 
(23) 



which are valid in the region 0.6 < T/T c < 0.99. Since the temperature of hadronic 
matter is generally larger than 0.11 GeV and smaller than T c , the meson masses and cross 
sections shown in the next paragraph in the temperature region are sufficient for studies 
concerned with hadronic matter. 

The wave function of quark-antiquark relative motion of the n and p mesons is taken 
to be the same as those of the other ground-state mesons. With m u = m d = 0.32 GeV 
obtained in Section 2, the experimental data of S-w&ve 7 = 2 elastic phase shifts for tttt 
scattering in vacuum [66-69] are reproduced with the potential in Eq. (17). Keeping 
quark masses at the values determined in the fit to the experimental mass splittings, 
the dependence of the transition amplitude A4a on temperature comes from the wave 
function, the potential in Eq. (17) and the meson masses. The transition amplitude 
contains also color, spin and flavor matrices that are not affected by the temperature. 
Since the meson masses depend on temperature, threshold energies for inelastic meson- 
meson scattering depend on temperature. The temperature dependence of the potential, 
the quark-antiquark wave function and the meson masses leads to temperature-dependent 
cross sections for inelastic meson-meson scattering. Unpolarized cross sections for the 
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seven reactions tttt -> pp for 7 = 2, KK -> 7T7T for 7 = 1, Kir* -> 7T7T for 7=1, 
vrTi" -> pTT* for 7 = 3/2, ttTT* -> p7i~* for 7 = 3/2, pTT -> pTT* for 7 = 3/2, and 
7r7\"* — > pK for 7 = 3/2 are plotted in Figs. 3-9. For convenient application the numerical 
cross sections are parametrized as 

, \ Cl r / i— r 

- V^o 



unpol 



'£o'" J 

(T'-'- = (l-i \ £-= j rxi) 



C2 



,'Vs- V*o\ 
■ <>■! \ £ ) ' >" !> 



C 2 (l 



(24) 



Time-consuming computations determine values of the parameters ai, 6i, Ci, a 2 , 62 and 
c 2 , which are shown in Tables 2 and 3. 

At the threshold energy of an endothermic reaction at a given temperature the mo- 
menta of final mesons in the center-of-mass frame equal zero. It is shown by Eqs. (12) 
and (13) that the absolute values of the momenta of initial and final mesons (| P \ and 
I P' |) increase as y^s increases. The rise of | P' \ causes a rapid increase of the cross 
section close to the threshold energy. The relative momentum p ab is a linear combination 
of P and P' . Thus, | p a b \ increases with increasing yfs, while ip a b{Pab) is reduced by the 
increase of | p a b \- The wave function ip a b(j>ab) and the absolute value of the transition 
amplitude | Ai& \ thus decrease with increasing y/s. This decreases the cross section with 
increasing yfs. The rising | P' \ and falling tp a b(p a b) produce a peak in the cross section 
near the threshold energy. At a higher temperature the constituents a and b have larger 
support of relative motion in coordinate space and so ipabipab) gets narrower in momentum 
space. This results in a cross section that decreases faster from the peak and forms a more 
narrow peak. The exception is that the width of the peak of a reaction at T/T c = 0.95 is 
equal to or slightly larger than the one of the same reaction at T/T c = 0.9. 

Cross sections for the seven reactions in Figs. 3-9 indicate that peak cross sections 
increase from T/T c = to 0.65 but decrease from T/T c = 0.9 to 0.95. While the tem- 
perature increases from zero, the potential plateau at large distances appears and moves 
lower. The Schrodinger equation with the potential produces increasing meson radii. 
Peak cross sections increase from T/T c = to 0.65 as the radii of initial mesons increase. 
While temperature goes to a higher value, bound states of mesons become looser due to 
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weaker confinement. At a temperature near T c the bound states are very loose. Even 
though mesons are easily broken in the reaction A(q 1 q 1 ) + B(q2q 2 ) — > C(q{q 2 ) + -D(?2?i) 
as T/T c — > 1, it is more difficult to combine final quarks and antiquarks into mesons 
through quark rearrangement. Hence, peak cross sections decrease from T/T c = 0.9 to 
0.95. On the one hand larger radii of initial mesons cause larger cross sections for the 
meson-meson nonresonant reactions, on the other hand looser bound states of final mesons 
lead to smaller cross sections. When temperature increases, both factors generate rising 
or falling peak cross sections. 

We may divide the seven reactions into three classes. The first class consists of the 
three reactions nn — > pp for 1 = 2, KK — > K*K* for I — 1 and nK — > pK* for 7 = 3/2 
which have spin for the two initial mesons and spin 1 for the two final mesons. The 
two maxima of the peak cross sections for the reactions are located at T/T c = 0.9 and 
0.75, respectively. The second class is just the reaction ttK* — > pK for 1 = 3/2 which has 
spin for an initial meson and a final meson and spin 1 for the others. The maximum 
of the peak cross section is located at T/T c = 0.85. The third class is comprised of the 
three reactions KK* -> K*K* for 1 = 1, nK* -> pK* for 1 = 3/2 and pK -> pK* for 
1 = 3/2 which have spin for an initial meson and spin 1 for the others. The maximum 
of the peak cross sections for these reactions is located at T/T c = 0.65 or 0.75. In the first 
or third class the meson-flavor dependence of masses distinguishes the peak values of the 
three reactions at a given temperature. In the first class the difference of the total mass 
of the initial mesons and the total mass of the final mesons is larger than the ones in the 
other classes. When T/T c changes from 0.6 to 0.99, the mass difference of the first class 
is reduced faster than the other two classes. The quicker reduction causes a more rapid 
increase of the peak value so that another maximum of the peak cross sections appears 
at T/T c = 0.9 in addition to the maximum at T/T c = 0.75. 

The potential in Eq. (5) results from the interpolation between the potential of per- 
turbative QCD at r < 0.01 fm and the lattice gauge results at large distances. The 
interpolation is in fact this procedure of adjusting the parameters A and E. We now 
begin to examine the sensitivity of mass splittings at T = GeV, temperature-dependent 
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meson masses and cross sections to the interpolation procedure by adopting two new sets 
of A and E. The first set (named Set I) is obtained by reducing A in Eq. (6) by 5% 
and increasing E in Eq. (7) by 10% while the second set (named Set II) is obtained by 
increasing A in Eq. (6) by 5% and reducing E in Eq. (7) by 10%. The two new sets give 
good fits to the lattice gauge results. Mass splittings at T = GeV and meson masses 
at five temperatures resulted from the two sets are listed in Tables 4 and 5, respectively. 
The change of the mass splittings or of the meson masses from Eqs. (6) and (7) to Set I 
or Set II is very small. The largest change in meson mass is 2.88% at T = 0.65T C . But 
such a change only leads to negligible changes in cross section. For example, the peak 
cross section obtained from Set I for KK — > K*K* for I = 1 at T/T c = 0.65 is 1.2735 
mb in comparison to the value 1.2686 mb obtained from Eqs. (6) and (7), and the peak 
cross section obtained from Set II for nK -> pK* for I = 3/2 at T/T c = 0.65 is 0.4787 
mb very close to the value 0.4930 mb resulted from Eqs. (6) and (7). Eventually, we 
understand that the mass splittings at T = GeV, the temperature-dependent masses 
and the cross sections are not sensitive to the interpolation procedure for the construction 
of the potential given in Eq. (5). 

In Figs. 3-9 we plot the average of the cross section in the prior form and the one in 
the post form. To see the uncertainty in this prescription, as an example, we draw in Fig. 
10 cross sections obtained in the prior form and in the post form for the reaction tck — > pp 
for 1 = 2. The solid curves in Fig. 3 are between the dashed curves corresponding to the 
prior form and the dotted curves corresponding to the post form. At a given temperature 
the dashed curve and the dotted curve almost overlap at the center-of-mass energies very 
close to or far away from the threshold energy. The post-prior discrepancy can be clearly 
marked by the difference of peak cross sections obtained in the two forms. We thus list the 
peak cross sections in Tables 6-8. The three tables are enough to display the discrepancy. 
Denote the peak cross section obtained in the prior (post) form by <jP™ r (crg^). To 
indicate the discrepancy, we define 

prior _ _post 

_ _ max w max (2^) 
X- prior , post \ / 

Cmax ~T~ C"max 
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Since o"P"° r > and crg^ > 0, the values of x are between -1 and 1. The smaller the 
absolute value of x, the smaller the discrepancy. If x is positive, cr^™ > crg^; otherwise, 
(jP™ r < o"max- The quantity x is al so presented in Tables 6-8. From Fig. 10 we note 
that the peak cross sections obtained in the two forms at a given temperature may not 
correspond to the same center-of-mass energy. 

5. Summary 

We have given the temperature-dependent central-spin-independent potential that in- 
terpolates between the perturbative-QCD potential with loop corrections at short dis- 
tances and the potential data offered by lattice gauge calculations at large distances. From 
the potential we obtain: (1) experimental mass splittings of the ground-state mesons with 
the same isospin when the masses of up, down and strange quarks are determined; (2) 
meson masses that decrease from T/T c = 0.6 to 0.99; (3) the wave function of quark- 
antiquark relative motion that occupies a larger volume at a higher temperature. In the 
quark-interchange mechanism we have obtained cross sections for seven nonresonant reac- 
tions 7T7T -> pp for 1 = 2, KK -> K*K* for I = 1, KK* -> K*K* for I = 1, nK ^ pK* 
for / = 3/2, irK* -> pK* for / = 3/2, pK -> pK* for / = 3/2 and ttK* -> pK for 
I = 3/2. The temperature dependence of the cross sections is determined by the tem- 
perature dependence of the potential, the quark- ant iquark wave function and the meson 
masses. Peak cross sections are affected by three factors: larger sizes of initial mesons at 
a higher temperature give larger peak cross sections, looser bound states of final mesons 
at a higher temperature lead to smaller peak cross sections, and a smaller total-mass 
difference of the initial mesons and the final mesons yields larger peak cross sections. The 
numerical cross sections are parametrized for future studies. 
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Table 1: Vacuum masses of 7r, p, K and K* and their spin- averaged masses. 





rn vp (GeV) 


m u (GeV) 


m p (GeV) 


rriKK* (GeV) 


m K (GeV) 


m K * (GeV) 


model 


0.7112 


0.2620 


0.8609 


0.8339 


0.5465 


0.9298 


experiment 


0.6124 


0.1396 


0.7700 


0.7929 


0.4957 


0.8920 
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Table 2: Values of parameters in the parametrization given in Eq. (24). 
reaction 

I = 2 7T7T — > pp 

I = 2 7T7T — > pp 

7 = 2 7T7T — > pp 

I = 2 7T7T — > pp 

7 = 2 7T7T — > pp 

7 = 1 7\~7T -> 7T7T 
7 = 1 7\~7T -> 7T*7r* 
7 = 1 TsTTT -> Tr*TT* 
7 = 1 7\~7T -> 7T*7\~* 
7 = 1 7\~7T -> 7T7T 
7 = 1 TOT* -> 7T*7\~* 
7 = 1 TsTTT* -> 7T*7\~* 
7 = 1 7\~7T* -> 7T*7\~* 
7 = 1 7\~7T* -> 7T*7\~* 
7 = 1 7\~7T* -> 7T7T 



T/T r 

1 C 


d\ (mb) 


6i (GeV) 


Ci 


do (mb) 


b 2 (GeV) 


Co 


0.65 


0.52 


0.13 


2.47 


0.07 


0.04 


0.42 


0.75 


0.86 


0.13 


3.11 


0.35 


0.04 


0.71 


0.85 


0.65 


0.0075 


0.38 


0.97 


0.12 


1.93 


0.9 


1.96 


0.009 


0.5 


0.43 


0.15 


1.58 


0.95 


1.03 


0.017 


0.47 


0.18 


0.28 


3.55 


0.65 


0.85 


0.11 


0.89 


0.4 


0.1 


0.44 


0.75 


1.42 


0.09 


1.14 


0.46 


0.15 


0.39 


0.85 


1 


0.06 


0.76 


0.47 


0.08 


0.22 


0.9 


1.3 


0.005 


0.63 


0.72 


0.04 


0.23 


0.95 


1.56 


0.006 


0.48 


0.3 


0.065 


0.31 


0.65 


1.71 


0.09 


0.58 


0.66 


0.1 


2.64 


0.75 


1.63 


0.09 


0.71 


0.83 


0.07 


1.61 


0.85 


0.88 


0.038 


0.57 


0.49 


0.1 


0.49 


0.9 


0.33 


0.006 


1.17 


0.57 


0.035 


0.21 


0.95 


0.51 


0.006 


0.47 


0.24 


0.05 


0.28 
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Table 3: Values of parameters in the parametrization given in Eq. (24). 



reaction 



T/T r 

/ C 


ai (mb) 


6i (GeV) 


Ci 


a? (mb) 


b 2 (GeV) 


Co 


0.65 


0.3 


0.01 


0.89 


0.48 


0.16 


1.45 


0.75 


0.44 


0.14 


1.43 


0.22 


0.1 


0.47 


0.85 


0.09 


0.006 


0.001 


0.41 


0.093 


0.96 


0.9 


0.76 


0.007 


0.65 


0.24 


0.05 


0.16 


0.95 


0.55 


0.012 


0.59 


0.09 


0.03 


0.15 


0.65 


0.72 


0.13 


1.29 


0.06 


0.06 


0.19 


0.75 


0.74 


0.11 


0.93 


0.05 


0.04 


0.42 


0.85 


0.22 


0.084 


1.95 


0.21 


0.08 


0.26 


0.9 


0.25 


0.007 


0.63 


0.19 


0.05 


0.24 


0.95 


0.16 


0.01 


0.56 


0.08 


0.07 


0.24 


0.65 


0.7 


0.001 


2.65 


1.1 


0.12 


1.2 


0.75 


0.69 


0.11 


0.81 


0.36 


0.09 


2.1 


0.85 


0.1 


0.11 


0.12 


0.36 


0.07 


0.84 


0.9 


0.27 


0.007 


0.53 


0.21 


0.05 


0.29 


0.95 


0.09 


0.045 


0.14 


0.09 


0.03 


4.54 


0.65 


1.68 


0.012 


0.53 


0.46 


0.06 


0.38 


0.75 


2.12 


0.01 


0.76 


0.66 


0.01 


0.09 


0.85 


1.82 


0.004 


0.48 


1.31 


0.01 


0.46 


0.9 


1.6 


0.005 


0.42 


0.29 


0.03 


0.68 


0.95 


0.58 


0.006 


0.4 


0.2 


0.04 


1.18 



J = I n K -> pK* 
I = | irK -> pK* 
I = | n K -> pK* 
I = | n K -> pK* 
I = | n K -> pK* 
I = | nK* -> pK* 
I = | nK* -> pK* 
1 = 1 nK* -> pis:* 
/ = | ttK* -> pK* 

1 = 1 7TK* -> piT* 

1 = 1 pK ->pK* 
I = \pK^pK* 
I = \pK -> pK* 
I = \pK ^ pK* 
I = \pK ^ pK* 
I = | ttK* -> pfsT 
J = | ttK* -> pfsT 

i = | vris:* _> 
i = I vris:* -> pis: 

i = I vriT* -> pis: 
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Table 4: Mass splittings in units of GeV at zero temperature. 

m p — m n rriK* — itik \m iJ j + |m<^ — m v 

Eqs. (6) and (7) 0.5989 0.3833 0.3622 

Set I 0.5890 0.3770 0.3530 

Set II 0.6090 0.3898 0.3666 
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Table 5: Masses of n, p, K and K* in units of GeV at various temperatures. 







T/T c 


m n 


m 


m K 


m K * 


Eqs. 


(6) and (7) 


0.65 


0.4042 


0.5776 


0.6161 


0.7271 


Eqs. 


(6) and (7) 


0.75 


0.3839 


0.4710 


0.5789 


0.6347 
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Table 6: Post-prior discrepancy 



reaction T/T c <g (mb) <g (mb) X 

I = 2 7T7T -> pp 0.5383 0.4443 0.0957 

7 = 2 7T7T -> pp 0.65 0.6208 0.5157 0.0925 

7 = 2 nn — > pp 0.75 1.2168 0.9064 0.1462 

I = 2nn^pp 0.85 1.2073 0.9010 0.1453 

7 = 2 7T7T -> pp 0.9 1.6551 2.2077 -0.1431 

7 = 2 7T7T -> pp 0.95 0.8296 1.2089 -0.1861 

7 = 1 KiT -> 7T*7T* 0.5386 0.7257 -0.1480 



I = 1 i^if -> 7T*7T* 0.65 1.2308 1.3510 -0.0466 

7 = 1 7T7T -> 7T*7T* 0.75 1.8543 1.8601 -0.0016 

I=IKK-*K*K* 0.85 1.7125 1.3213 0.1289 

I = 1 for -> 7T7T 0.9 1.5737 2.2874 -0.1848 

7 = 1 7T7T -> K*K* 0.95 1.5242 2.0874 -0.1559 

7 = 1 TTTT* -> 7T7T 0.9475 0.9425 0.0026 

7 = 1 TTTsT* -> 7T7T 0.65 2.1129 2.6930 -0.1207 

7 = 1 TTTT* -> 7T*7s:* 0.75 2.3603 2.6151 -0.0512 

I = 1KK*^K*K* 0.85 1.4742 1.1657 0.1169 

7 = 1 TT7<* -> TTTsT* 0.9 1.6670 0.0377 0.9558 

7 = 1 TTTC* -> 7sT*7T* 0.95 1.3724 0.0163 0.9765 
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Table 7: Post-prior discrepancy 



reaction T/T c <^ r (mb) <g (mb) ^_ 

I — \ ttK — > pK* 0.3430 0.3615 -0.0263 

I = \ixK^pK* 0.65 0.5111 0.4830 0.0283 

I = \ixK ^ pK* 0.75 0.7084 0.6268 0.0611 

I = \ixK^pK* 0.85 0.6149 0.4901 0.1129 

I = \nK-+pK* 0.9 0.7691 1.1911 -0.2153 

I = I ttK ^ pK* 0.95 0.4953 0.7650 -0.2140 

I = \nK*^pK* 0.5415 0.4310 0.1136 

J = | TriT* -> pfT* 0.65 0.7573 0.8279 -0.0445 

1 = 1 tiK* -> pif* 0.75 0.8269 0.8191 0.0047 

I = \^K*^pK* 0.85 0.5147 0.4022 0.1227 

J = | ttK* -> pK* 0.9 0.6862 0.1273 0.6870 

I = \ rrK* -> pfT* 0.95 0.3985 0.0823 0.6576 

I = \ pK ^ pK* 0.6318 0.5553 0.0644 

I = \pK^pK* 0.65 1.1125 1.1724 -0.0262 

I = \ pK pK* 0.75 1.1026 1.0686 0.0157 

I = \pK~^pK* 0.85 0.5808 0.4368 0.1415 

I = \ pK ^ pK* 0.9 0.8565 0.0061 0.9859 

I = \ pK pK* 0.95 0.4294 0.0080 0.9634 
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Table 8: Post-prior discrepancy 
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Figure 1: Temperature-dependent potential. From top to bottom the temperatures corre- 
sponding to the different data sets [25] are T/T c = 0.58, 0.66, 0.74, 0.84, 0.9, 0.94, 0.97. 
Solid curves stand for the parametrization fitted to the data at these temperatures. 
= 2.154 fm" 1 = 0.425 GeV. 
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Figure 3: txtx — > pp cross sections for / = 2 at various temperatures. 
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Figure 4: KK — > K*K* cross sections for / = 1 at various temperatures. 
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Figure 5: KK* — > K*K* cross sections for / = 1 at various temperatures. 
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Figure 6: irK — > pK* cross sections for I — 3/2 at various temperatures. 
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Figure 7: nK* — > pK* cross sections for I = 3/2 at various temperatures. 
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Figure 8: pK — > pK* cross sections for I — 3/2 at various temperatures. 
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Figure 10: irn — > pp cross sections for 7 = 2 in the prior form (dashed curves) and in 
post form (dotted curves) at various temperatures. 
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